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Abstract. We report the first time-resolved photometric and spectroscopic optical observations of the X-ray source 
RX J2 133.7+5 107, identified in the ROSAT survey. A clear persistent optical light pulsation is discovered with fast photometry 
at a period of P w =(570. 823+0.013) s which we associate with the spin period of an accreting white dwarf. Radial velocity curves 
of the strong emission lines show modulation with a period of Pn=(7. 193+0.016) hr, identified as the orbital period. These ob- 
servations establish that the source is a member of the intermediate polar class (IPs) of magnetic cataclysmic variables. With 
only 4 IPs with longer orbital periods, RX J2 133. 7+5 107 is among the widest systems. It is a unique IP with an orbital period in 
the middle of the so-called (6-10)hr IP gap and it shows a significant degree of asynchronism with a ratio P w /Pn of 0.02. When 
attributed to the motion of the white dwarf, the emission lines orbital modulation yields a mass function of f„, = (1.05+0.21) 
xlO -2 M Q which, for a probable inclination i < 45° and a white dwarf mass M w a = (0.6 - 1.0)M Q , corresponds to a secondary 
mass M s > (0.27 - 0.37) M e . 
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the WD rotational period. In an effort to elucidate the nature of 
The X-ray source RX J2133.7+5107 was discovered from the RX J2133.7+5107, we obtained optical observations to search 
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k*" ■ 1. Introduction of channelled accretion in these systems is a strong pulsation at 

CS ; ROSAT Galactic Plane Survey and identified as a relatively for a signature of the WD rotation and its orbital motion 
bright (m B ~ 16) cataclysmic variable (CV) (Motch et al. 1998, 
hereafter M98). Its X-ray hardness ratios indicate a hard X- 

ray source, similar to the intermediate polar (IPs) class among 2. Photometric Observations 

cataclysmic variables (CVs). IPs are Magnetic Cataclysmic 

Variables (MCVs), where a white dwarf (WD) with a rela- RXJ2133. 7+5107 was observed from August 23fh to 26th, 

tively strong (B < 5-10MG) magnetic field accretes material 2003 at the 1.5m Loiano telescope (Bologna, Italy) equipped 

from a late-type secondary companion overflowing its Roche- with the three channel photometer (TTCP) with simultaneous 

lobe. The accretion flow is channelled towards the polar re- acquisition of target, comparison star, and sky. Observations 

gions where a strong shock develops above the WD surface (see were done without filter with a photometer efficiency peaking 

Patterson (1994) for a review), which means that the signature around 3800 A and a spectral range of -3000-6000 A (Silvotti 

et al. 2000). The integration time was set to 30 s in all ob- 

~ ; ~ : T ,, , , servations. Sky conditions were photometric with a seeing of 

Send offprint requests to: J.M. Bonnet-Bidaud, email: bonnet- r^„s 
bidaud@ceafr during all nights, except on August 24 (seeing of 2 ) 

* LJTon observations obtained at the Haute-Provence when observations were also interrupted due to bad sky con- 
Observatory (France) and at the Loiano Observatory (Italy), op- ditions for -1.3 hr. A diaphragm of 17" was used for all three 
erated by the Istituto Nazionale di Astrofisica ' channels. RX J2133.7+5107 was estimated at B~16mag us- 
** Ta bles 2 and 3 and Figures 4-5- 6 are only available in electronic ing a nearby comparison star. The log of the observations is 
form at http://www.edpsciences.org reported in Tabled The photometric data were reduced using 

Draft 1.0, February 5, 2008 



J.M. Bonnet-Bidaud et al.: RX J2133.7+5 107: a new long period Intermediate Polar 



- Aug. 23, 2003 



■ Aug- 24. 2003 



- .iS::,. . i 



^ Aug. 26, 2003 



A 1 * V.. 



I 



3.4 3.5 

Time (HJD 2452875.0 + ) 



Fig. 1. RXJ2133. 7+5107 photometric light curves in relative 
white light amplitudes with a 30 s resolution. A ~ 10 min. pulsa- 
tion is clearly seen on all three nights, superposed on an overall 
long term variation. 
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Fig. 2. Power spectrum of the photometric data (pre-whitened 
from the low frequency trend) around the spin fundamental 
(bottom) and first harmonic (top) frequencies. Note the differ- 
ent power scale for the two ranges of frequencies. Left panels: 
direct Fourier transform of the time series, central panels: syn- 
thetic power spectrum using the best fit periods (see text). Right 
panels: time series "clean" power spectra after removing the ef- 
fect of the data window. The location of the fundamental (of) 
and orbital sideband (a> - Q.) are shown by vertical dotted lines. 
For the first harmonics, the small peak around ~ 302 day 1 is a 
residual 1-day alias that has been imperfectly removed. 



standard procedures including sky subtraction, extinction cor- 
rection, and relative photometry. The resulting light curve is 
shown in Fig^ 

The photometric time series were analyzed performing a 
Discrete Fourier Transform (DFT) of the relative amplitudes, 
after detrending from the low frequency underlying variations. 
A significant power was present in two groups of frequencies 
inside the overall range (140-160) day -1 (see Fig. |2j. This 
structure closely resembles the one observed in the optical 
power spectra of many IPs with the presence of the spin fre- 
quency (oj) and its orbital sideband (10 - Q), where £2 is the 
orbital frequency. The same structure was also found, though 
with much lower intensity, at twice these frequencies, which 
are identified as the first harmonics (see Fig. |2] top panels). 
Given the strong influence of the data window, a number of 
aliases were also introduced that prevent any clear identifi- 
cation of the true peaks. To remove this effect, we used the 
CLEAN algorithm ( Robert s"et al. 19871 . and results are also 
shown in Fig.|2](right panels). The most significant peaks were 
found at co — 151.3+0.1 day -1 with its orbital sideband at u-Q. 
= 148.0+0.1 day" 1 and harmonics at 2u = 302.6+0.1 day" 1 , 
where the error bars were determined from the peaks'half 
widths. The fundamental and the first harmonic of the spin 
match strictly, while for the beat frequency, the weak peak ob- 
served at 296.6+0.1 day" 1 , is slightly offset from the expected 
beat harmonic (2a> - 2Q) frequency. This discrepancy is most 
likely due to the very small amplitude (full amplitude 0.8%) 
and imperfect removal of the data window effect. 

As the presence of a sideband beat period may significantly 
affect the period determination, the detrended data were fitted 
with a four frequency sinusoidal function that includes the fun- 



damental and first harmonic of the spin and beat periods. The 
best periods were then found at : P w = 570.8227+0.013 s and 
P<y-a= 583.768+0.02 s. This fit was used to derive the best spin 
ephemerisforRXJ2133.7+5107 as: 

T™ = 2 452876.94785(3) + 0.0066067(2) E (1) 

where T™ ax is the heliocentric arrival time of the sine curve 
maximum. A quadratic fit on the (O-C) residuals does not yield 
a significant period derivative with an upper limit P < 2.0 10~ 6 
s/s. 

We checked the consistency by producing a synthetic DFT 
of the composite sinusoidal function, including the four most 
significant frequencies above. Results are shown in Fig[2](cen- 
tral panels). The good agreement between the observed and 
synthetic spectrum confirms that the spin and beat modulations 
have been determined accurately. The data, pre-whitened from 
the low frequency trend, were folded at the derived spin pe- 
riod and the pulse shape is shown in Fig [3] The pulse has a 
full amplitude of ~ 8% and a quasi-sinusoidal shape with a less 
steep rise than the decay, indicative of the presence of the first 
harmonic. 

From the photometric data, the measured beat frequency 
yields an orbital frequency of Q = (3.357±0.006)day~ 1 corre- 
sponding to Pn = 7. 150+0.015 hr, in agreement with the orbital 
period found from spectroscopy (see below). The limited num- 
ber of nights and the probable contamination from a likely non- 
periodic variability do not allow a more precise determination 
of the orbital period. 
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Table 1. Log of photometric and spectroscopic observations. 
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Obs ervjitiori s 


Range 


Res. 


Date 




UT(start) 


Exposure (min) 


Photometry 


White light 


30 s 


2003 Aug 


23 


21:15 


362 




White light 


30 s 


2003 Aug 


24 


20:09 


167 




White light 


30 s 


2003 Aug 


25 


00:17 


172 




VV 111 IL llgllL 


JU a 








HO 7 


Spectroscopy 


3600-7200 


5.7 A 


1998 Aug 


24 


20:35 


165 




3600-7200 


5.7A 


1998 Aug 


26 


20:53 


210 




4300-6000 


2.8 A 


2004 Jul. 


20 


20:39 


270 




4300-6000 


2.8 A 


2004 Jul. 


21 


21:08 


270 




4300-6000 


2.8 A 


2004 Jul. 


22 


20:34 


300 



RXJ2133 + 51 




5000 5500 
•velenglh ( A) 



Rotational F 



Fig.4. The mean optical spectrum of RX J2133.7+5107 ob- 
tained in July 2004, showing typical strong Balmer and Helium 
emission lines. Note the significant CIII-NIII and the presence 
Fig. 3. Pulse shape of the 57 1 s spin modulation. A significant of diffufje ^5^^ bands (DIB) and Na absorption, blended 
distortion from a pure sine shape is visible due to the presence w j tn ^ jj e j (5875A) emission 
of a first harmonic. 



3. Spectroscopic observations 

RXJ2133. 7+5107 was observed at the Observatoire de Haute 
Provence (France) at two epochs in August 1998 and in July 
2004 (see Table 0. Long slit spectra were obtained with the 
Carelec spectrograph (Lemaitre et al. 1990) attached to the 
Cassegrain focus of the 193 cm telescope and using a CCD 
EEV (2048x1024 pixels) detector of 13.5/im pixel size. In 
1998, the 133 A/mm grating was used, with a slit width of 2" , 
leading to a wavelength coverage of 3600-7200 A at an FWHM 
resolution of ~ 5.7 A. In 2004, spectra were obtained with the 
67 A/mm grating and a slit of 2" width, covering the range 
4300-6000 A at a resolution of 2.8 A. Exposure times were ei- 
ther 15 min or 30 min in 1998, and 30 min in 2004. In 1998, 
clouds were present on the last night (26 Aug). For both nights, 
the seeing was greater (3-5 ) than the slit (2" ). In 2004, ob- 
servations were obtained in nearly photometric conditions with 
only thin cirrus during the last night (22 Jul). Seeing was typi- 
cally of 2.5" to 3" . 

Standard reduction was performed in the ESO-MIDAS 
package, including cosmic removal, bias subtraction, flat-field 
correction, and wavelength calibration. The wavelength cali- 
bration was checked on sky lines that were found within 1.6 A 
and 0.7 A from their expected wavelengths in 1998 and 2004, 
respectively. All radial velocity measurements were corrected 



from the Earth motion and from the small instrumental shifts 
measured on the OI 5577 A line and times were converted in 
the heliocentric system. Flux calibration was performed using 
the standards BD+28 421 1 and BD+33 2642. 

The mean optical spectrum of RX J2133.7+5107 (Fig.|4) 
is typical of magnetic CVs with strong emission lines of 
the Balmer series, Hell (4686A and 5411 A) and CIII-NIII 
(4655A), superposed on a relatively blue continuum. In 2004, 
when the flux calibration was reliable, the mean B flux was es- 
timated at ~ 3 x 10~ 15 ergs cm _2 s _1 A , which corresponds to 
a B magnitude of ~ 15.8, similar to what was reported in the 
identification paper by M98. The characteristics of the main 
emission lines are given in Table 13 Similar equivalent widths 
were observed in the two observations in 1998 and 2004, but 
the intensities of the lines change by a factor 2. 

The radial velocities were measured for the main 
lines. They showed a significant modulation of ~ (30-70) 
km s _1 amplitudes that were analysed to constrain the orbital 
period. The exposure of the individual spectra (30 min.) were 
close to a multiple of the spin period so that no strong influence 
from the spin modulation was expected, but the period search 
was significantly affected by the 1-day sampling of the observa- 
tion. The best period was determined from an x 1 periodogram 
where the radial velocity measurements were folded at trial pe- 
riods and tested against a sine modulation. The most accurate 
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Table 2. Equivalent widths, FWHM and fluxes of the strongest lines 



Line 


EW(A) 


FWHM(A) 


Flux(*) 


EW(A) 


FWHM(A) 


Flux(*) 






August 1998 






July 2004 




Htf 


2.9(3) 


7.2 


56 




. 




H y 


3.9(2) 


7.5 


73 


4.2(2) 


7.3 


136 


HeI4471 


1.0(1) 


7.5 


17 


0.9(1) 


6.2 


30 


CUJ-NIU(+) 


2.0(1) 


21.3 


36 


2.3(2) 


26.9 


71 


Hell 4686(+) 


5.9(1) 


10.6 


110 


5.7(2) 


10.5 


182 




6.1(3) 


8.5 


109 


6.4(3) 


8.0 


201 


Hell 5411 


1.5(3) 


12.8 


23 


1.4(2) 


15.2 


42 


Hel 5875(+) 


2(2) 


8.1 


32 


2.0(2) 


7.9 


55 


H, 


16.5(4) 


10.7 


252 









Errors bars on last digit in parentheses. 

FWHM: corrected for instrumental resolution (typical errors of 0.6A in 1998 and 0.3A in 2004) 
* Line flux in units of ergs cirT 2 s -1 are indicative because of marginal photometric conditions 
+ Measurements derived from Gaussian fits. 



determination was obtained for the strongest line (H^ ) in the 
longest observation of 2004 (Fig|5j. The minimum^- 2 corre- 
sponds to a period of Pn=7. 193+0. 016hr, where the error bar 
is at a 99% level, computed for two independent parameters. 
A secondary minimum, also observed at ~ 10.3 hr, is a 1-day 
alias that can be excluded at a level better than 3. 10~ 4 from an 
F-test. A similar value was obtained independently in 1998 by 
making use of the velocities of the strongest H a line, but the 
accuracy on the period value does not allow us to keep track of 
the phase between the two observations, so the orbital period 
cannot be refined. The orbital ephemeris for RX J2133. 7+5107 
is determined as : 

T n = 2453208.1428(28) + 0.29973(67) E (2) 

where Tn is the predicted heliocentric time of the blue-to- 
red radial velocity transition. The radial velocities of the main 
lines, folded at the best orbital ephemeris above, are shown in 
Fig-El an d the parameters of the corresponding best sine fits are 
given in Table[3] Data are from the 2004 observations obtained 
with the highest spectral resolution, except for H a , which was 
only observed in 1998. For this last line , the orbital phase was 
arbitrarily shifted, as the ephemeris were not accurate enough 
to link the 1998 data. The Balmer lines showed similar ~ 
70 km s velocity amplitudes, while the Hell lines displayed 
lower value with no significant phase shift. Inspection of the 
line intensities, EW, and widths folded at the orbital period do 
not reveal significant modulation, except for a very small vari- 
ation in width for H r and . 

4. Discussion 

The optical characteristics of RXJ2133. 7+5107 clearly point 
toward an IP nature of the source. Photometric data show a 
clear, persistent 571s period, which we interpret as the spin 
period of the accreting white dwarf. The presence of a spin 
first harmonic and the detected beat periodicity consistent with 
a ~ 7.2 hr orbital period are features common to many IPs 
systems, signaling a probable two-pole accretion with signif- 
icant reprocessing inside the binary system (see lWarner 19951 



12000 




6 8 10 12 

Period (hr) 

Fig. 5. x 1 periodogram of the radial velocities of in 2004 
(black line) and H a in 1998 (grey line). The same best value 
is found in the two observations for Pn ~ 7.2 hr. Secondary 
minima are 1-day (2004) and 2-day (1998) aliases 



Table 3. The parameters of the orbital modulation of the emis- 
sion lines' radial velocities. 



Lines 


y ( kms ') 


K(kms _1 ) 


Phase (*) 


H y 


+ 12.7 (1.4) 


66.5 (2.0) 


0.003 (0.030) 




+ 10.6 (0.4) 


67.4 (0.6) 


0.000 (0.009) 


HeII( 4686 A) 


- 14.7 (0.9) 


42.2(1.2) 


0.024 (0.029) 


HellX 5411 A) 


- 34.2 (5.2) 


33.5 (7.4) 


0.021 (0.220) 


H, 


- 22.8 (0.4) 


74.8 (0.6) 


0.0 (fixed) 



* Phase of the blue-to-red zero crossing 



1-cr error bars into parentheses 



|Buckley 2000[ >. The orbital period is confirmed independently 
by modulation of the emission lines, which yields a more 
precise value of P = (7.193+0.016)hr. This is among the 
longest orbital periods found in intermediate polars as, out of 
~40 systems, only four confirmed IPs (AEAqr, V1062Tau, 
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Fig. 6. Radial velocities of the main lines folded with the orbital 
period P n =7.193hr. 



RXJ1548-4528 and GKPer) known so far have wider orbits 
JNorton et al. 2004> . 

With a longer period than 7 hr, a secondary of a rela- 
tively early type (up to K0) can be expected. Such stars usu- 
ally show a significant sodium line Na D (5889-5896A) in 
absorption. This line is indeed detected in the spectrum of 
RXJ2133. 7+5107, partly blended with the Hel (5875A) line, 
at the resolution used in our observations (see Fig.0}. However 
the Na D line observed in RX J2 13 3. 7 +5 107 is at least partly 
interstellar, indicative of both significant distance and absorp- 
tion. This is confirmed by the presence in the spectrum of the 
diffuse bands at DIB(/14430) and DIB(/t5780) similar to what is 
observed in some other CVs (see |Thorstensen & Taylor 200 1) . 
The equivalent widths of the (^4430) and (/15780) DIBs are 
measured at 0.7A and 0.3A respectively, and are both con- 
sistent with a reddenning of Eb-v ~ 0.3 dHerbig l975| 
ISnow et al. 20021 . If purely interstellar, the Na Dl (/15890) 
component would have indicated a much higher reddenning 
of Eb-v ~ 1-1 (Munari & Zwitte r~1997l l. inconsistent with the 
above value and indicating that the secondary also contributes 
to the line. From the measured equivalent widths, the star con- 
tribution can be estimated to be less than 50 %. 

The Na D equivalent widths were measured on representa- 
tive stars from the Jacoby et al. (1984) catalogue with spectral 
types from K0 to M0 and compared to the EW upper limit mea- 
sured in RXJ2133. 7+5107. If the secondary contributes for 
less than half of the observed Na D, the contribution of a K0, 
K5, or M0 star to the optical flux is respectively less than 39%, 
37%, and 9%, which translates into a lower limit to their dis- 
tance of 1200, 650, and 600 pc respectively, assuming standard 



absolute magnitudes (Allen 1970). X-ray observations of the 
source by ROSAT do not give better indications of absorption 
and distance. From the observed ROSAT hardness ratios and 
the measured value HR1 =0.90+0.06 (M98), using the HR1-N H 
distribution for CVs in Fig. 4 (M98) gives a rough estimation of 
Nh ~ 1 x 10 21 cirT 2 with, however, a large uncertainty. The 
measured ROSAT flux and the count conversion given in M98 
yield a [0.1-2.4 keV] luminosity of L x ~ 1.5 x 10 32 (d/500pc) 2 
ergs s _1 . 

A complementary constraint on the system can also be ob- 
tained from the velocity amplitude of the emission lines, if 
they are assumed to be associated with the white dwarf. The 
mean velocity of the Balmer lines yields a mass function of 
f m = (1.05+0.21) xl0~ 2 M Q . For an inclination of i < 45°, 
suggested by the absence of double-peaked lines, and a white 
dwarf mass M w( j = (0.6 - l.O)M , , the secondary mass is 
M s > (0.27 - 0.37)M o , consistent with a star at least more 
massive than M2V. 

The relatively long orbital period of RX J2133.7+5107 
makes it fall into what has been called the "IP gap" 
JSchenker et al. 20041 . a lack of IPs into the (6-10) hr orbital 
period interval. From the list of 39 IPs in No rton et al. 2 004 
all are either shorter than 6.14 hr or longer than 9.87 hr, this 
last value being the one newly determined for 1RXS J 1548- 
4528 by de Martino et al. (2005). Though it could still be due 
to an observational bias, this may indicate a particular state of 
evolution. 

With a spin-to-orbit period ratio of 0.02, RXJ2133. 7+5107 
is also a relatively highly asynchronous system with only 8 
IPs showing lower values. This suggests a young system still 
far from equilibrium. For magnetic CVs, Norton et al. (2004) 
have computed the equilibrium rate when the spin-up accretion 
torque is balanced by the braking effect of the magnetic torque, 
using a detailed plasma-magnetic field interaction. For a value 
of Po,/Pfi= 0.02 and an orbital period of 7.2hr, the equilibrium 
is reached for a magnetic moment of fi ~ 10 33 G.cm 3 , here as- 
suming a mass ratio q = 0.5, but the variation is less than 10% 
for q in the range 0.3< q <0.7. The magnetic moment is prob- 
ably much lower in RX J2 13 3. 7 +5 107 , since accretion will 
only take place if the Alven radius (r m ) is lower than the coro- 
tation radius (r co ) (see lWarnerT9^5l l. For RXJ2133. 7+5107, 
r co = (0.87 - 1.03) x 10 10 cm (for M wd in the range (0.6-1.0) 
M ). For// ~ 10 33 G.cm 3 , the (r m < r co ) condition will require 
an accretion rate that is higher than (2.8-3.5) xlO 17 g.s , so 
much too high for the observed X-ray luminosity. Assuming 
a more canonical value of M ~ 10 16 g.s~', the magnetic mo- 
ment then has to be lower than (1.2-1.9) xlO 32 G.cm 3 ; there- 
fore RXJ2133. 7+5107 is most probably an IP with a weak 
magnetic field. A definite confirmation of its IP nature will be 
possible with the detection of strong X-ray pulses expected at 
the 571 s period. 
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